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Abstract. The perturbations of the gravitational field due to the mass distribution 
of an absolute gravimeter have been studied. The so called Self Attraction Effect (SAE) 
is crucial for the measurement accuracy, especially for the International Comparisons, 
and for the uncertainty budget evaluation. Three instruments have been analysed: 
MPG-2, FG5-238 and IMPG-02. The SAE has been calculated using a numerical 
method based on FEM simulation. The observed effect has been treated as an 
additional vertical gravity gradient. The correction (SAC) to be applied to the 
computed g value has been associated with the specific height level, where the 
measurement result is typically reported. The magnitude of the obtained corrections 
is of order 1 x 10~^ m/s^. 
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1. Introduction 

Modern transportable Absolute Gravimeters (AGs) measure the local value of the free- 
fall acceleration g by reconstructing the trajectory of a falling object in vacuum. The 
mass of the parts constituting each AG apparatus (such as laser interferometer, vibration 
isolation system, vacuum chamber, etc.) are sources of an additional gravitational field, 
which can systematically perturb the motion of the flying object. The so called Self 
Attraction Effect (SAE) could be non negligible because the actual absolute g values 
are estimated with uncertainties of the order of few yuGal (1 yuGal = 10^ ms~^) and 
masses of few kilograms placed at few centimeters from the test body can generate 
fluctuations of the gravitational field of several microgals. 

To calculate the SAE, a detailed study on three different gravimeters MPG-2 [Ij, 
FG5-238 ^ and IMGC-02 [3] has been performed, using a Finite Element Method 
(FEM) simulation to calculate the contributions of each part of the gravimeters. Then 
the Self-Attraction Correction (SAC) to be applied to the g value has been calculated 
for each AG, taking into account the peculiar model used to compute the free-fall 
acceleration from the reconstructed free-fall trajectory. 

The knowledge of the SAE and the relative SAC is crucial for each measurement 
carried out by the AGs, because it can improve the accuracy as well its uncertainty. 
It becomes very important when the measurements are used to calculate the Key 
Comparison Reference Value (not physically known) during the dedicated International 
Comparisons. 

In the following section a brief description of the three instruments is presented. 
After that, the adopted method to evaluate the SAE and the results obtained for the 
AGs are illustrated. In the last part, the SAC values are calculated separately for the 
three AGs. 

2. Absolute gravimeters 

Three different transportable AGs have been studied: 

• the MPG-2, designed in Germany by the Max Planck Institute for the Science of 
Light, prototype instrument; 

• the FG5-238, actually developed and commercialised by the U.S.A. Micro-g LaCoste 
Incorporation (the results coming from this AG can be considered for to the other 
instruments of the same series, i.e. FG5-2a;a;); 

• the IMGC-02, developed in Italy by the Istituto Nazionale di Ricerca Metrologica 
(INRIM), prototype instrument. 

For each of them, the measurement of the g value is obtained by the reconstruction 
of the trajectory of a corner-cube prism (or a hollow retroreflector) which moves 
vertically in vacuum. The IMGC-02 takes into account for both the rise and fall motions 
of the flying object, while the other two instruments measure the acceleration during 
the free-fall motion only. 
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Automated systems are employed to centre, launch and receive the object event 
by event with nominal rates of about 0.02-h0.05 Hz during data taking sessions of 
several hours. An interferometric system is implemented to measure time and distance 
coordinates of the trajectory using a visible laser beam. The interferometer measures the 
distance between the corner-cube flying prism and a second one mounted on the quasi- 
inertial mass of a vibration isolation system used to protect it from ground vibrations. 
A detailed description of the tree AGs can be found in [H |2l |3]. 

The distribution of the heaviest and nearest parts of each instrument and the path 
of the flying object have been considered to characterise the perturbation gravitational 
field T{Z). Each gravimeter can be essentially divided in three main parts: 

• read-out electronic case: it can be easily moved farther away (up to about 1 m) 
from flying object trajectory, so the SAE from this source is negligible; 

• measuring system: supporting tripod, seismometer or super-spring system with its 
support, several detectors and the interferometer, in this case the SAE is mainly 
due to the masses of each single parts; 

• launch system: the vacuum chamber with the dropping mechanism and its basis 
whith all the accessories, in this case the effect must to be analysed in details 
because those objects are very close to the flying object trajectory. 

The geometry of the last two parts have been drawn using the COMSOL® 3D 
module and it is shown for the three AGs in figure [H Screws, cables, holes are not 
simulated because their influence is negligible for the SAE. 

3. Self-Attraction Effect 

To evaluate the gravitational field perturbation of each single part of an AG, an accurate 
knowledge of the geometry and the mass distribution of the source is needed. Due to 
the complexity of the single parts (edges, different materials, non regular shape), a FEM 
simulation has been preferred instead of a mathematical modelling. 

3.1. Numerical method 

The COMSOL® software [1] has been adopted. A package devoted to gravitational 
effects is not implemented in the original software. The electrostatic module is then 
used, exploiting the analogy between gravitational and electrical interactions already 
used for this or similar purpose in [5l [6] . Defining a conversion parameter as 



with G the universal gravitational constant and Eq the electric permittivity constant in 
the vacuum, the module of the gravitational field is given by 



K = -AttGeo 



(1) 




(2) 
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Figure 1. Geometry of the MPG-2 (top), FG5-238 (middle) and IMGC-02 (bottom) 
absolute gravimeters, drawn using the COMSOL® 3D module. The two main systems 
are distinguishable: measuring (coloured) and launch one (transparent). Only parts 
interesting for the gravitational effect analysis have been simulated. 
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where pe and pm are respectively the electrostatic and mass density, V the volume of 
the part. 

The component along the direction orthogonal to the floor on which a gravimeter 
is located (called Z coordinate) of the gravitational field (Tz) has been considered to 
evaluate possible effects on the measurement of the g value. The other two components 
along the X and Y axes can influence the flying object trajectory introducing rotation 
or shift components. This effort has not been treated in this study for two principal 
motivations: i) the effects are related to the mechanical features of each AG, as the 
corner-cube centring or the asymmetry in the launch system; ii) the perturbations are 
expected to be negligible with a respect to the one along the Z coordinate. 

For these reasons, the origin of the reference frame has been set on the ground floor, 
with the Z axis directed upwards and crossing the centre of the flying object. Modern 
AGs measure the g value at the height values between about 0.4 and 1.3 m, so the Tz 
value is calculated along a straight line from 0(0, 0, 0) to A{0, 0, 1.5 m). 

The FEM simulation consists of three main steps: geometry draw, mesh 
implementation and equations solution. 

In the first step, the parts of the AG are drawn and their mass values implemented. 
Two different approaches can be used for each part: 

• if its shape is regular (such as cylinder or parallelepiped) and the mass distribution 
is almost homogeneous, the mass density and volume real values are introduced in 
the software; 

• in case of complex shapes, e.g. different parts with screws and holes, or 
inhomogeneous mass distribution, as a seismometer, it is approximated to a regular 
solid with equivalent volume and mass is simulated. 

In this study, seismometers, frames, launch chambers, vacuum pumps have been 
simulated using the second approach, the other components have been simulated with 
the real shape. 

The parts so defined are then embedded in a empty volume where the field can be 
propagated. The shape is spherical with the radius value about one order of magnitude 
larger than the sizes of the AGs, in order to have edges far from the centre of the 
instrument, so as to minimise the boundary effects [6]. 



In the second step, the mesh geometry and size are implemented. Such parameters 
must be tuned to have the maximum resolution of the gravitational field in the range 
of interest. 

To validate the FEM simulation parameters, i.e. the mesh size and the boundary 
sphere, a study of two simple objects has been preliminary performed. A steel sphere 
centred at 0(0,0,0) with radius ri = 0.12 m has been located together with an 
aluminium sphere centred at P(0, 0, 0.6 m) and r2 = 0.1 m. The radius of the boundary 
sphere is Tq = 3 m. In figure |2] the relative difference between the values of the 
gravitational field (only component along Z coordinate at X=0 and Y=0) obtained from 
the theoretical model and the FEM simulation is plotted. No increment of the residuals 
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can be observed far from the masses {Z ^ A), indicating that no effect due to the 
boundary conditions is present. The fluctuations are caused by the mesh discretization. 
The residual values are distributed around zero with the root mean square below 2%, 
as shown in the inset plot. 

In the last step of the method, the equations of the problem are numerically solved. 
In this case the only parameters can be set are referred to the CPU consumption for 
the algorithm, with modifications of the results below the 0.01%. 




0.4 0.6 
z coordinate [m] 



Figure 2. Comparison between gravitational field values (component on Z coordinate, 
X=Y=0) obtained by using the theoretical formulas and the FEM simulation with a 
so called finer mesh. 



The main contribution to the uncertainty of the simulation are due to the 
approximations made in the first step. In order to estimate it, several simulations 
of different complex parts (e.g. seismometer, launch chamber) have been performed 
varying the ratio between the sub-part masse^. Using such approach, the standard 
uncertainty of the SAE computed using the FEM method has been estimated to be 
MsAE = 0.1 yuGal. 



3.2. SAE results 

The method described above has been applied to the three AGs with equal FEM 
parameters. A mesh called finer, which is automatically adapted to the edges and 
the contact points, has been implemented. The gravitational field has been calculated 

J In other words, an object of known mass is constituted by sub-parts of unknown mass. Simulations 
with a ratio between the sub-parts set to 1/3, 1/3, 1/3 is performed. Then it is repeated with different 
ratio values as 1/8, 3/8, 1/2 etc. 
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in the whole boundary sphere (ro = 3 m). The component Z of the gravitational field 
along the line centred with the flying object trajectory has been extracted for the global 
simulation. The contributions due to the measuring and the launch systems have been 
separately estimated. 

In flgure [3] the global and the single results are shown for the three AGs. 

For the MPG-2, the largest values are located at Z ^ 0.8 m and Z ~ 1.2 m, due 
to the presence of the launch chamber basis. In the range of the free-falling path, an 
average SAE of about 1.5 /iGal is estimated and the behaviour appears monotone. 

Also for the FG5-238 the main effect can be observed at Z ^ 0.8 m, but in this 
case the free-fall range is larger than the previous instrument 0.2 m). The SAE in 
such range is below 2 /iGal with an average value of 1.3 yuGal. 

For the IMGC-02, two bumps can be observed at Z ^ and Z ~ 0.9, due to the 
presence of the launch chamber basis and the seismometer with its support respectively. 
In the range of rise and fall path, the mean influence of the SAE is below 0.2 /iGal 
(absolute value). 

4. Self- Attraction Correction 

The computed above SAE is nearly linear in the measurement range of every AG 
(flgure [3]). Such linear (over height) gravity variation corresponds to the constant self- 
attraction gravity gradient 7sae with a some constant offset zsae at a chosen origin. The 
superposition property of gravitational flelds from several sources allows to treat the SAE 
in an AG similar and in addition to the conventional vertical gradient 7 of the Earth 
gravity fleld. Then the observed height-varied free-fall acceleration is approximately 
given by 



The z-axis in ([3]) is directed toward the centre of Earth (opposite to the Z axis in 
the previous flgures), and the signs of go and 7 are chosen to be positive downwards 
(a freely-faUing object accelerates downwards, magnitude of the acceleration increases 
towards the Earth surface). This section considers methods to convert the computed 
SAE into the SAC for the particular AG. 

4-1- Effect of the gravity gradient 7 

The linear free-fall motion model for the unperturbed free-fall in a homogeneous gravity 
fleld reads [8] 



where zq and vq are the initial position and velocity at t = 0, and g^ is a value of the 
free-fall acceleration at the level z = 0. When the registered free-fall motion trajectory 



g{z) = {go + -fz) + (iJsAE + 7sae2;)- 



(3) 



z{t) = zo + vot + -got^, 



(4) 
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Figure 3. Result of the FEM simulation for the absolute gravimeter MPG-2 (top), 
FG5-238 (middle) and IMGC-02 (bottom). Self- Attraction Effect versus Z coordinate. 
Total and single contributes of launch and measuring systems are shown. The 
approximate operating range of every instrument is identified by the rectangular in 
the figure and it is enlarged in the inset plot. 
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Zi{ti) is constructed from more than three time-distance coordinate pairs, the standard 
least-squares estimate on qq in a model (jl]) is given by the functional |9] 



g = G(ti, Zi) 



YjtiZi 



N Uti T.t1 
T.ti Stf Stf 
St? St? Stf 



(5) 



N Uti 
T,ti St2 
St? St^ 

where sums comprise all N time-distance coordinates ti and Zi. 

In a case of the non-homogeneous gravity field characterized by the constant gravity 
gradient 7, the approximated free-fall motion model reads 



t 



t' 



z{t) = zAl + -i- \ + t;o U + 1 + ^?o ( - + 7 



t' 



(6) 



^2 ' '24^ 

The least-squares approximation of the model (jlj) instead of ([6]) gives a biased 
estimate g, which equals to an instantaneous gravity value at the specific point below 
the level ^ = 0, called the effective measurement height h^s [S [lOl [TTl [12] : 

fl- = fl-o + 7^cff- (7) 

The remarkable property of the /leff is that the solution g (equation ([5])) at this point 
does not depend on the constant gravity gradient and its uncertainty. It can be found 
as [9l[l2] 



ieff 



^0 +^^0 



G{ti,ti) G{ti,tf) 

—^ + 9^^, (8) 

where the parameters estimated zq, vq and g are obtained from the least-squares 
approximation of the model (jl]), and the functional ([5]) is modified by replacing Zi 
with if at Vq or with at g. An asymptotic (A^ ^ 1) expression^ for the hes can be 
found in P[ini[n]. 

When the fitted model ([6]) contains the gravity gradient 7 as an external known 
quantitjjjj], the obtained least-squares estimate go is automatically translated from the 
hes to the level 2; = as 

9o = 9-lKs- (9) 

Since the considered absolute gravimeters employ different models, it results in 
different methods of computing the SAC. 



4.2. SACforMPG-2 

The MPG-2 is a multiple-level free-fall absolute gravimeter with the linear free-fall 
motion model (jll) [I], therefore, the measurement result is reported at the hew- A 
straight-line approximation of the SAE (figure HJ left) is given by 

^(^) = ^SAE + 7SAE (^ - /ieff)- (10) 

§ The effective measurement height is usually counted from the start [zi) of the measurement distance 
interval H [lOl [12]. To make it compatible with the h^s ([8]) used in this report, one should add a 
distance interval from the level z — (e.g., optical centre of the resting retro- reflector) to Zi. 
II When the model ([6]) does not contain the gravity gradient at the term zq |2], the reference level 
corresponds to z = zq instead of z = 0. 
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The parameters zsae (intersect) and 7sae (slope) can be found by the least-squares 
approximation of the computed SAE with the associated standard deviations Mzs^e and 
^7sAE- ^ range of approximation should match the actual measurement distance interval 
H = {zf — Zi), used to compute ([5]). 

Top of drop (z = 0) 




Figure 4. Self-attraction effect over a measurement range H of tlie MPG-2, FG5-238 
and IMGC-02. Open grey dots for the SAE computed in previous section; red lines 
for the approximations by the first order polynomial. Solid squares indicate a value of 
the SAC (opposite sign) at the corresponding reference height. Computed corrections 
are shown in table [1] 

According to (IHl) and ([7]), the measurement result at the hes is given by 



g = 90 + Zsae + (7 + 7sae) Ks- (11) 

Since this estimate does not depend on the constant gravity gradient 7tot = 7+7sae5 
the correction is given by 

= -ZSAE, (12) 
with a standard uncertainty of 

UAg=U-isAE- (13) 



The the computed SAC is of -(1.6 ± 0.004) /iGal. It should be added to the 
result reported at the level of i^inst — ^eff- 

4.3. SAC for FG5-238 

The FG5-238 is a multiple-level free-fall absolute gravimeter with the linear free-fall 
motion model ([6]) (without the gradient term at ^o) P], therefore, the measurement result 
is obtained at the level zq (with further reduction of (—7^0) it is usually transferred 
to the rest position of the falling object, so called "top of drop" [2J). Combining (fTTIl 
and (IHl), the measurement result is given by 

g = go + Zsae + 7sae hes. (14) 

It has to be corrected by 

Ag = -(^sAE + 7sAE Ks)- (15) 
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The standard uncertainty of this correction is given by: 

UAg = U^sAE + ^ISAE heS- (16) 

Note that the standard uncertainty ( fT6l) under the same conditions is larger 
compared to the uncertainty (fT3ll . This is a particular outcome of a general disadvantage 
to include a known gravity gradient into a model ([6]) and to report the measurement 
result at the top of drop instead of reporting it at the effective measurement height. 

If a value of the effective measurement height is not available, the computed 
SAE can be approximated by the straight line with the origin placed to the starlit 
of the measurement interval H. Then similar to (fT5|) . a correction can be obtained by 
extrapolation to the reference level (e.g. top of drop): AgsAE = — (0.7± 0.004) /iGal. 
It should be added to the result reported at the instrumental reference height i^mst 
(figure m centre). Alternatively, the SAC at the level of -ffinst — ^eff is given by 
-(1.1 ±0.003) /iGal. 



4.4. SAC for IMGC-02 

The IMGC-02 is a multiple-level rise-and-fall absolute gravimeter with the non-linear 
model, which contains the gravity gradient 7 as an unknown parameter [13] : 



z{t) = Zo + ga 



{t-t,f {t-t,f , {t-t,Y 



V 7 + 7 



(17) 



6 ' 24 

where ta is the time at the apex of the trajectory parabola and ip the friction coefficient 
of the residual air in the launch chamber. 

Since this model is non-linear in the parameters ga and 7, the superposition principle 
for the perturbed parameter estimates is not valid, and a posteriori correction for 
a time (or height) varying perturbation is not applicable in a strict sense. Instead, 
a straightforward approach to modify the registered time-distance coordinates before 
the non-linear least-squares adjustment might be used by taking into account a 
priori computed SAE. However, obtained results might not be extendible to other 
measurement conditions or changed parameters of an instrument. Then such a 
modification is necessary at every new measurement session. 

In the IMGC-02 a measurement result is usually reported at the best measurement 
height Zh below the apex of the parabola [13]: 



gb = ga + izb, (18) 

where ga and 7 are directly estimated from the trajectory data using the model (1T71) . 
and Zf) is computed from the elements of the variance-covariance matrix and also depends 
on a quantity of the upper portion of the trajectory removed from the least-squares fit. 
With typical parameters of the IMGC-02 the best measurement height is close the 



^ Shifting an origin to the border of the range of approximation gives rise to the standard error of the 
extrapolated parameter, which is, however, neghgibly small for the particular curvature of the SAE. 
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Table 1. The self- attraction corrections for three absolute j 


gravimeters. 




MPG-2 


FG5-238 


lMGC-02 


refer, height above floor [m] 


1.154 


1.292 


0.475 


measur. range H = {zf — Zi) [m] 


0.094 


0.192 


2x0.174 


effective height /iofi[ni] 


0.043 


0.080 


0.029 


start velocity v{zi) [m/s] 


0.31 


0.35 


1.84 


SAC (linear approx.) [/^Gal] 


-1.6 ± 0.1 


-0.7 ± 0.1 


+0.6 ± 0.1 


SAC (mean value) [/iGal] 


-1.5 ± 0.3 


-1.3 ± 0.3 


+0.2 ± 0.4 



effective measurement height /icfr = H/Q derived for the rise-and-fall AG with a hnear 
free-fall motion model and without removing the upper portion of the trajectory [9]: 

'-'-Is <^«' 

with H being a total distance interval traveled by a falling object twice. 

To evaluate a magnitude of the SAE on the measurement result, we assume that 
at the level the estimate gb (HlSi) is not correlated with the estimate 7 |13]. Then 
the bias is given only by the constant term of the SAE, derived from the least-squares 
adjustment of the model (ITOl) with an origin chosen at the level Zb- Consequently, the 
SAC can be evaluated by f|T2l) with the standard uncertainty f|T3|) . The SAC is then 
computed as +(0.6 ± 0.0001) /iCal. It should be added to the result reported at the 
level of i^inst - Zb (figure IH right). 

Table [1] contains the computed corrections for the considered AGs. The standard 
uncertainty from the least-squares approximation is upgraded to 0.1 yuGal, as derived 
from the FEM simulation. In addition, a crude mean- value estimate on the SAC is also 
considered (corresponding to the average SAE values with an opposite sign). It might 
be seen, that this simple estimate is a satisfactory measure of the SAC in a case when 
the known gravity gradient is not included in the free-fall motion model, or when the 
measurement result is reported at the effective (or the best) measurement height. 

Reporting the g measurement result at the top of drop should not be associated 
with the averaged SAE. A proper SAC requires further extrapolation outside the actual 
measurement range H (figure HI middle). 



5. Conclusions 



Three Absolute Gravimeters (AGs) have been studied in order to calculate the Self- 
Attraction Effect (SAE) due to the masses of their single parts, using a FEM simulation 
method. The correction (SAC) to be applied to the final measure of g has then been 
calculated for each AG. It should be pointed out that besides the observed SAE, the SAC 
generally depends on the method of a free-fall acceleration measurement, adopted free- 
fall motion model, number and method of data location and portion of the reconstructed 
free-fall motion trajectory used in the least-squares adjustment. 
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The MPG-2 presents the largest SAE with an average value of (+1.5 ± 0.1) /iGal 
and a SAC can be applied to the reported g value. For the FG5-238 an averaged SAE 
of about (+1.3 ±0.1) /iGal has been found, consistent with the results estimated by [7]. 
Also in this case a SAC of (—0.7 ± 0.1) /iGal can be applied to the measurement result 
reported at the top of drop. For the IMGC-02 the averaged effect (—0.2 ± 0.1) /xGal 
results to be negligible with a respect to the combined standard uncertainty. However, 
a a priori approach to implement the SAC has been discussed. 

Such results can be related to the measurements performed during all the 
previous and future comparisons, especially the International Comparisons of Absolute 
Gravimeters (ICAGs). The most recent published results of an ICAG (i.e. 2005 [H]) 
have been chosen as a reference in the following. 

In the case of the MPG-2 (not present at ICAG 2005) the computed SAC exceeds the 
preliminary declared uncertainty of 0.5 /xGal due to this effect. Therefore, in the results 
and uncertainty budgets of MPG-2, reported at previous comparisons, it is reasonable 
to replace this contribution of 0.5 /iGal with the obtained SAC. 

About the IMGC-02, the application of the SAC should not significantly change 
the g values with a respect to the combined standard uncertainty of about 4.3 /iGal, 
which already included a contribution of 0.3 /iGal for SAE |14j . 

The results found for the FG5-238 can be considered for the other AGs of the same 
series (i.e. FG5-2a;x) which represent the most part of the instruments participating 
to the comparisons. Usually the g measurements performed using FG5-xxx did not 
present a SAC and only few participants considered the SAE in the uncertainty budget 
evaluation. At ICAG 2005 [14|, a contribution of 0.1 /iGal has been estimated for the 
FG5-108 and the FG5-213 and it is has been used to calculate the so called conventional 
uncertainty for all the FG5-a;a;a; [15]. Hence, a revision of the values of the previous 
comparison should be useful, also because the abundance of FG5 instruments strongly 
influences the Key Comparison Reference Value. 

The presented approach to evaluate the SAE and the SAC is applied to three AGs, 
but it has general validity. For this reason, the whole procedure can be applied to other 
absolute gravimeters, knowing its the peculiar features as geometry, mass values and fit 
model. 
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